Spontaneous Initiation of Detonations in White Dwarf 
Environments: Determination of Critical Sizes 



Ivo R. Seitenzahl 1 ' 2 , Casey A. Meakin 2,3 ' 4,5 , Dean M. Townsley 2,4 , Don Q. Lamb 3,4 , James 

W. Truran 2 - 3 ' 4 ' 6 - 7 



ABSTRACT 

Some explosion models for Type la supernovae (SNe la), such as the gravita- 
tionally confined detonation (GCD) or the double detonation sub-Chandrasekhar 
(DDSC) models, rely on the spontaneous initiation of a detonation in the de- 
generate 12 C/ 16 material of a white dwarf. The length scales pertinent to 
the initiation of the detonation are notoriously unresolved in multi-dimensional 
stellar simulations, prompting the use of results of ID simulations at higher 
resolution, such as the ones performed for this work, as guidelines for deciding 
whether or not conditions reached in the higher dimensional full star simulations 
successfully would lead to the onset of a detonation. Spontaneous initiation relies 
on the existence of a suitable gradient in self-ignition (induction) times of the 
fuel, which we set up with a spatially localized non-uniformity of temperature 
- a hot spot. We determine the critical (smallest) sizes of such hot spots that 
still marginally result in a detonation in white dwarf matter by integrating the 
reactive Euler equations with the hydrodynamics code flash. We quantify 
the dependences of the critical sizes of such hot spots on composition, back- 
ground temperature, peak temperature, geometry, and functional form of the 
temperature disturbance, many of which were hitherto largely unexplored in the 
literature. We discuss the implications of our results in the context of modeling 
of SNe la. 
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Introduction 



Type la supernovae (SNe la) are believed to be thermonuclear explosions of accret- 
ing white dwarfs, powered by the energy liberated in the fusion of the initial composition 
to more tightly bound nuclear species, often all the way to nuclear statistical equilibrium 
(NSE). The explosion mechanism of SNe la is still unknown, and the large variation of 
their peak luminosity and spectral properties allows for a range of explosion mechanisms. 
Pure deflagration models have a lot of unprocessed carbon and oxygen at low velocities and 
generall y turbulently mixed la ye rs of nuclear materi al burned to differing degrees of com- 
pletion ( Reinecke et al. 2002al lbl; iGamezo et al. 2003 ). This is more or less in contradiction 
to observations, which show a layered st ructure of the ejecta, and a lack of low velocity 
low and intermediate mass elements 
remaining viable exp 



(e.g. 



Khokhlov et al. 1997 



onation (GCD) fjPlewa et al. 



Meakin et al. 
1990l : 



Branch! 120041 ) . Therefore, many of the currently 
qsion models, such a s deflagration to deto nation transition (DDT) (e.g. 



Gamezo et a 



Fink et al. 



2004 



2005 



Townsley et al. 



2007 



Ro pke et al.ll2007ah. gravitat i onally confined det- 



Plewa 2007; Jordan et al. 



2008), or double detonation sub-Chandrasekhar (DDSC) (ILivne &: Glasner 



2008 



20071 ) have detonations as a necessary ingredient. 



The spatial scales relevant for the initiation of these detonations are unresolved in even 
the most computationally expensive multi-dimensional sim ulations, and recen t major efforts 



i n the modeling of the SNe l a exp losion mechanism li ke iFink et al.l (120071 ). iRopke et al. 



( 12007bl ) , iTownsley et al.l (120071 ) , and IJordan et al.l (120081 ) all rather similarly decide whether 
or not a detonation is launched based on the peak temperature reached above a certain 
density in a computational cell. To make the call in favor of or against a successful detonation, 
one commonly checks whether or not a ID reactive hydrodynamics calculation starting with 
a hot spot consisting of a linear temperature profile, with peak temperature equal to the 
hottest temperature that obtains in the simulation, leads to a detonation at that density for 
reasonably small v alues of the radius of the te mperature inhomogeneity, as is done in the 
pioneering work of iNiemeyer fc Woosleyl (119971 ). 



Such detonation conditions derived from ID unresolved simulations, which assume a lin- 
ear temperature gradient, a cold isothermal ambient medium, spherical geometry, constant 
density and fuel concentration, and an initially stagnant velocity flow field, are clearly an 
oversimplyfication of the initiation problem. Therefore, currently employed procedures of 
calling the outcome of a supernova simulation based on external detonation conditions (in- 
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cluding the ones presented here) should be taken with a grain of salt and prudence should be 
exercised in their application. In this paper, we illustrate the uncertainties in the detonation 
conditions by taking a closer look at the consequences of relaxing the stringent assumptions 
about the temperature profile and geometry. 

Since existence of a supersonic induction time gradient due to a temperature and/or 
fuel concentration gradient (and to a lesser degree density) is the essential ingredient in the 
formation of the detonation, it is crucial that the high temperatures be reached in unburned 
fuel and that the gradients be not too steep. For unresolved temperature gradients, the spa- 
tial variation of the fuel concentration, which is an important ingredient in the determination 
of the induction time gradient, is unknown. Even when assuming the best case scenario of 
100% fuel concentration, one should allow for the possibility of different functional forms 
of the unresolved temperature gradient, the effects of which we explore in this paper. As 
we demonstrate in this work, for the same density, peak and background temperatures a 
gaussian temperature profile may not lead to a detonation where a linear profile does, even 
though more energy was placed into the hot spot region. Similarly, for some choices of en- 
vironmental parameters a gaussian temperature profile leads to detonation where a linear 
profile does not, changing the conclusions about the success or failure to detonate of some 
simulations. Additionally, as we will show, the temperature of the medium surrounding the 
hot spot (the background temperature) can also have a large effect on the outcome and 
should be taken into account as an additional parameter. This work furthermore includes 
a study of the effects of composition (different mass fractions of 16 0, 12 C and 4 He) on the 
critical sizes. 

In section 2, we give an overview about detonations, their spontaneous initiation, and 
dependences on environmental parameters. In section 3 we describe our method for deter- 
mining critical radii for linear temperature profiles and critical decay constants for different 
exponential temperature profiles. Section 4 contains a presentation of the results and trends 
with varying environmental parameters. Section 5 closes with the conclusions. 



2. Theory of detonations 
2.1. Detonation primer 



In spite of the easy availability o f good introductory and review material on detonations 
( iFickett fc David Il979l ; IClavinl 120041 . i.e.), for the benefit of the reader we shall give a brief 
overview here. A detonation in its simplest form is a shock that advances supersonically into 
a reactive medium, behind which chemical or nuclear reactions proceed. At least part of the 
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region where exothermic reactions take place must be in sonic contact with the shock to inject 
energy and prevent the deto nation from fail i ng through dis sipation. The Chapman- Jouguet 
(CJ) model for detonations (IChapmanlll899l ; |Jouguetlll905l ) is a one-dimensional model that 
describes the detonation in the limit of infinite reaction rate. In the state immediately 
behind the shock, the reactions are assumed to have progressed to completion. The Rankine- 
Hugoniot jump conditions, which express the conservation of mass, momentum, and energy 
flux across the shock front are simply modified by including the compositional change and 
the associated energy source term across the shock in the conservation equations. Assuming 
the upstream material to be at rest, the conservation equations accross the detonation front 
are: 



PiD 



p 2 (D - u 2 

(ei + dq) + — + \d 2 = s 2 + — + l(D 



P1+P1D' 



U 2 



U 2 



(1) 

(2) 
(3) 



Pi 2 p 2 2 

where D is the detonation speed, dq is the specific energy liberated in the burn, and the other 
quantities take on their usual meanings. Using the specific volume V{ = 1/pi, and eliminating 
the fluid velocity of the ash state, u 2 , from equations [1] and [21 we get an expression for the 
so called Rayleigh0 line: 



Vi-Vt 



(4) 



Using eq. H] and either CD or [2] one can easily eliminate D 2 and (D — u 2 ) 2 from the energy 
conservation equation to get an expression for the Hugoniot curve: 

2 ( £l + dq - e 2 ) = (P 2 + P X ){V 2 - V x ) (5) 

Neumann! (Il942l ) and iDoerin improved the one dimensional 



ZePdovichl (11940 ) 



von 



description of detonations in what has become known as the ZND structure. They model a 
detonation as a one dimensional leading shock wave moving at the detonation speed trailed 
by a reaction zone of finite width in which energy is released and the fuel transforms into the 
burning products. The flow is assumed steady and planar. Solving the Rankine-Hugoniot 
equations without allowing for compositional change across the shock gives the thermody- 
namic conditions immediately behind the shock, the so called von Neumann state. From 
the von Neumann state on downstream the spatial variations of the hydrodynamic and 
thermodynamic variables in the reaction zone are determined by differentially following the 
compositional transmutations and the associated heat release. The ZND model has the ad- 
vantage that it gives a prescription for calculating the width of the detonation wave (i.e. the 



Also known as the Mikhel'son (sometimes spelled Michelson) line after the Russian physicist. 
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distance between the shock and the final state of the ash). For one particular detonation 
speed D C j, the Rayleigh line is tangent (i.e. one point of intersection) to the Hugoniot curve 
in the P — V plane, and the detonation is of so called CJ-type. In this case, the reactions 
terminate at the sonic point, that is, D — u 2 = c s , where c s is the sound speed in the ash. 
If the final ash composition is known, this uniquely determines the propagation speed. It 
turns out that this so called CJ-speed corresponds to the smallest possible detonation speed, 
and many detonations fall into this category (but see also section [276]) . For a detonation 
speed D > D C j, the Rayleigh line intersects the Hugoniot curve in two distinct places. The 
two points of intersection correspond to the "weak" or "under-compressed" and the "strong" 
solution. Strong detonations, which display the higher pressure increase and compression of 
the two solutions, have the reaction products moving subsonically with respect to the shock. 
Such detonations are oft en also referred to as "overdriven" , with overdrive factor / 
(e.g. 



Hwang et al J 12000( 1 . 



- D • 
Do J ■ 



It is known that no steadily propagating weak detonation is possible (e.g. lFickett Davis 



19791 ). Supersonic and shockless weak detonations however do occur as a transient during 
the early stages of the initiation of detonations from shallow induction time gradients where 
they play an integral part in the formation of the ultimately steadily p r opagating strong o r 
CJ detonation structure fcapila fc Rovtburdl[l~989l : Ishort fc Doldlbopj iKapila et al.ll2002h . 



While many of the observed properties of detonations can be explained by the simple one 
dimensional CJ and the slightly more advanced ZND models of steady state detonations,it is 
clear today that real detonations are very complex multi-dimensional structures that are at 
most steady state in an average sense. The structure and instabilities of fully developed det- 
onations are very complex and hav e received considerable attention from researchers in both 



the astrophysical community (e.g. iGamezo et al- 



and e specially the combustion community (e.g. 



19971 : IShort fc Stewart! Il998l : IShort & Sharpd 120031 ) . 



19991: ISharpelll999l: iTimmes et alihoOObh 



Erpenbeck 



1964 



1966 



197 



0; 



Clavin et al. 



2.2. Initiation of detonations 

The question of how to form or initiate a detonation is similarly still an active area of 
both experimental, theoretical, and numerical research. The different ways of detonation 
formation can be grouped into two distinct categories: 



1. Direct initiation 

2. Spontaneous (initially shock- less) initiation. 



- 6- 



Direct initiation involves a blast wave or an otherwise formed shock propagating into a reac- 
tive medium at a speed exceeding the CJ speed for the conditions in the fuel. The resulting 
overdriven shock reaction- zone co mplex trans itions or relaxes under the right conditions into 



a self supporting detonation (e.g. lBodyill997l ). 



Spontaneous initiation, which this work focuses on, involves shock formation due to a 
spatial gradient in the initial conditions of the induction times of the fuel. It should be 
noted that spontaneous initiation may still require an external shock to "precondition" the 
fuel such that auto-ignition may proceed from then on. The difference to direct initiation 
is that it is not the original shock that transitions into the leading shock of the detonation, 
but rather a new shock that forms when the nuclear fuel runs away coherently. The gradien t 
mechanism of the initiation of a detonation was first proposed by IZel'dovich et al.l (ll970l ). 
The principal feature of this mechanism is the presence of a gradient in induction times 
(self-ignition delay times) that leads to a supersonic reaction wave. In the spontaneous wave 
picture, the phase velocity of the burning front is given by: 



Vph 



t dr i s 



(6) 



This wave transitions to a detonation when its phase velocity becomes equal to the velocity 
of a Chapman- Jouguet (CJ) detonation in that material. Zel'dovich's spontaneous wave con- 
cept, which ignores non-linear gas dynamical evolution and derives the burning wave speed 
from the initial conditions of the initiating center, describes the initiation of a detonation 
very well for shallow gradients. For steep gradients near cri ticality, however, t he nonlinear 
gasdynamical effects are large and essential to the problem (IKapila et al.ll2002l ). 



Lee et al.l (119781 ) build on Zel'dovich's idea and publicized the g eneralized picture of 



shock wave amplification through coherent energy release (SWACER). iBartenev fc Gelfand 
( 120001 ) summarize the key features of the mechanism as follows: 



Gas layer(s) with minimal induction time(s) ignite(s) first 

Shock propagates to adjacent layer(s), which are on threshold of ignition. 

Shock initiates instantaneous explosion of layer, which strengthens shock. 



In the astrophysical context o f thermonuclear supernovae. lBlinnikov fc Khokhlovl (1l986t ) 
and iBlinnikov fc Khokhlovl (119871 ) were among the first to discuss the initi ation process of 



Khokhlov et al. 



deton ations from a temperature gradient in degenerate 12 C/ 16 matter 
dl997h studied induction time gradient initiations for detonations that were not solely due 
to temperature but also fuel concentration. Pioneering work for determination of critical 
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radii in the context of SNe la ha s been perfo r med by lArnett &: Livnei (11994 ) , followed by 
Niemeyer fc Woosleyl (119971 ) and iRopke et al.l (l2007bl ). This paper extends their work to 



include a study of the dependences of the critical gradients on composition, background 
temperature, peak temperature, geometry, and functional form of the temperature distur- 
bance. Next, we discuss each of these briefly to motivate the components of our study. 



2.3. Fuel density 

The higher the fuel density, the smaller are the characteristic length scales for deto- 
nation. This is essentially due to the fact that thermonuclear reaction rates scale as the 
product of the number densities of the reactants (i.e. r it j = riiUj < av >ij). As a result 
of the higher reaction rates, burning time scales are shorter and the width of a fully de- 
velo ped self-support ed 12 C/ 16 detonation wave therefore decreases with increasing density 
(e.g. 



<ed sen-support ed u/ u detonation wave tnereiore decreases witn increasing density 
Khokhlovi[l989l l The critical length scales (e.g. critical radii) for the initiation process 



are correspondingly smaller as well. At low densities, the critical radius becomes larger than 
the density scale height of the star, leaving constant density simulations rather meaningless. 
Nevertheless, the density (for a given choice of peak and background temperature) at which 
the critical radius becomes equal to the size of the white dwarf can be considered a very 
conservative estimate of the lowest density at which a gradient initiated detonation could be 
successfully launched. A fiducial density of 10 7 g cm -3 is adopted here for most cases, but 
densities as low as 10 6 g cm -3 and as high as 3 x 10 7 g cm -3 are considered. 



2.4. Temperature profiles 

Past work in the determination of critical sizes of gradient initiated detonations in WD 
environments all considered linear temperature gradients and a cold background. Typically, 
for a given spatial extent of the temperature disturbance, the smallest peak temperature 
that would still lead to detonation in a cold surrounding medium was determined. A slightly 
different approach is taken here: For a given peak temperature the smallest size of the heated 
region that still leads to detonation is determined (see section [3]). 

The choice for the linear profile of the temperature gradient seems arbitrary. For un- 
resolved simulations, a hypothetical Gaussian or exponential temperature gradient on the 
sub-grid scale is at least as plausible as a linear one. Whether or not there is considerable 
dependence of the detonation conditions on the choice of functional form for the temperature 
gradient is therefore a question well worth investigating. Since the initiation site is possibly 
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embedded in an environment where even the asymptotic background temperature has been 
heated to an appreciable temperature (by a shock for example, or by compression), a range 
of ambient background temperatures is examined here. 

The initial temperature profiles of the simulations consist of a hotspot with peak tem- 
perature T max falling off to an ambient temperature To (see fig. [1]). The functional form of 
the temperature perturbation is either linear 



T max ~ Tmax R To r for r < R, 



T(r) 



for r > R, 



(7) 



exponential 



T(r) = {T rnax - T ) exp{-r/R) + T 



Gaussian 



T(r) = (T max - T ) exp[-(r/R) 2 ] + T , 



(9) 



higher order Gaussian like (hereafter referred to as glO) 



T(r) = (T max - T ) exp[-(r/i?) 10 ] + T ( 



(10) 



or a top hat 

T max for r < R, 
T(r)={ (11) 
T for r > R. 

Depending on whether the problem is setup in planar or spherical geometry, r is the distance 
from a reflecting boundary or the origin respectively. 



2.5. Composition 



SNe la progenitors are likely 12 C/ 16 white dwarfs. The 12 C to 16 ratio in a massive 
white dwarf at the end of core 4 He burning is not very well constrained. The large uncertainty 
of one of the most important nuclear reaction rates, 12 C(a,7) 16 0, translates into a range of 
possible 12 C to 16 ratios. In this work, we explore the sensit ivity of the critica l sizes on the 
12 C to 16 ratio for a range of allowed ratios (see fig. 4 from iDommguez et al.l (120011 )). 
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The details of the accretion process towards a Chandrasekhar mass white dwarf in the 
single degenerate channel are not solved, but the possibility of a significant layer or at least 
admixture of 4 He in the outer layers of the massive white dwarf remains. We would like to 
stress that the rate of the 12 C(a,7) 16 reaction proceeds much faster than carbon burning. 
Therefore, adding 4 He into the mixture will lead to a more reactive medium and significantly 
decrease the critical radii. In this work, we explore the sensitivity of the critical sizes on the 
possible admixture of 4 He. 



2.6. Geometry and curvature 



Planar detonations in 12 C/ 16 Q mat t er are known to be of pathological t ype for densities 



2 x 10 7 g cm" 3 (IKhokhlov! Il989l : ISharpel Il99ft iGamezo et all ll999T ). Endothermic 



P > 

photo-disintegration reactions give rise to a frozen sonic point in the reaction zone of the 
detonation at which the local sound speed equals the velocity of the flow in a frame in which 
the leading shock of the detonation is stationary. Material downstream of the sonic point 
is moving supersonically away form the shock front and is out of causal contact with the 
leading edge of the detonation. The sonic point, also known as the pathological point, lies 
for CO WDs generally in the oxygen or silicon burning layer. The detonation speed of a 



pathological detonation is an eigenvalue of the steady equations (e.g. lYao &: Stewartll2006l ) 



and it is larger than the detonation speed of the corresponding CJ detonation. 

Steady planar detonations in 12 C/ 16 O at lower density are of the C J- type, frozen subsonic 
throughout all the burning zones with the sonic point located at the end of the reaction zone 
where nuclear statistical equilibrium is reached. It is well known, however, that even a 
small amount of curvature can significantly influence the propagation speed and structure 
of detonations in SNe la environments. The curvature causes 12 C/ 16 detonations to be of 
pathological type even for densities as low as p = 5 x 10 6 g cm -3 . For curvatures larger than 
1.5 x 10~ n cm -1 , the pathological sonic point falls near the end of oxygen burning, greatly 
affecting the length scales of carbon, oxygen, and silicon burning stages of the detonation 
( )Sharpell200ll ). The curvature constraint means in practice that all detonations that initiate 



from a localized hot spot near the surface of a white dwarf for densities p ~ 10 7 g cm 3 are 
of pathological type. 



He fc Clavinl (119941 ) show for direct initiation, that non-linear curvature effects cause the 



critical radius R c , in this case defined as the radius correspondi ng to the sm allest possible 
igniter energy at which the shock wave of the Sedov problem (jSedovl 119591 ) transitions or 
relaxes to a wave propagating at the CJ-velocity of the mixture, is about 10 times larger for 
spherical geometry when compared to planar. It is intuitively clear that the critical radii 
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for spontanous initiation in ID spherical geometry will be larger compared to ID planar 
geometry. To quantify the differences critical radii in both geometries are determined for the 
same composition, density, background temperature and peak temperature. 



2.7. Multi-Dimensional Effects 



Real detonations are not one dimensional objects. They rather exhibit a multi-dimensional 
cellular structure with complicated internal substructures that includes transverse shocks, 
triple points (points of ma ximum pressure where transverse shocks inter sect the Mach stems) 
and weak incident shocks (ITimmes et al.ll2000bl ; iFickett fc David Il979l ). The cellular struc- 
ture behind the detonation front results in pockets of incompletely burned fuel which leads 
to a reduction in ene rgy input compared t o the one dimensional case and results in a reduced 
detonation velocity (IBoisseau et al.lll9961 ). 



While the cellular structure is important for the width of a self sustaining steady deto- 
nation and the resulting pockets of inhomoge neously burned fuel will likely leave a signature 
in the spectra of SNe la (IGamezo et al.lll999l ). it is unclear how it affects the critical radii for 
initiation of detonations via the gradient mechanism. An attempt to address this question 
was made by determining critical radii for a few select cases in 2D. Significant differences 
in the critical radii when going from ID to 2D were not noted. The results however, should 
be considered inconclusive, since the simulations started from idealized perfectly smooth 
initial conditions. Without initial perturbations, the time scale for the cellular instabilities 
to develop from numerical round off errors is large compared to the time scale for oxygen 
ignition. It is conceivable that slightly more realistic noisy initial conditions will facilitate 
faster growth of the cellular instability and change the picture altogether. Furthermore, it 
is conceivable that one has to spatially resolve the carbon burning layer to see a substantial 
manifestation of higher dimensional effects. This question certainly deserves further detailed 
study, which unfortunately goes beyond the scope of this article. 



Critical Size Determination 



Past efforts in determining minimal cri tical length scal e s for t he initiation of deto- 



nations in w 
and 



(1997) 



lite dwarf environme nts include lArnett &: Livnd (119941 ) , iNiemeyer fc Woosley 
Ropke et al.l (j2007bl ). All of these works are one dimensional simulations 



Arnett &: Livnel ^ 1994 ) determined critical radii in 12 C/ 16 material by determining the 



smallest spheres with a peak temperature T n 



2.0 x 10 9 K and a background tempera- 
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ture To = 2.0 x 10 s K and a top hat temperature profile that would lead to a runaway. A 
successful detonation was declared if the heated material was completely burned within 10 
sound crossing times (1000 time steps). Although a thermonuclear runaway of the heated 
region is a necessary ingredient for the formation of a detonation, it is by no means sufficient. 
Additionally to the runaway, a shock-reaction zone complex needs to form and survive. The 
functional form of the temperature profile and the extent of the heated region both play a 
significant role here. 



To illustrate, a hot spot in ID planar geometry with a top hat temperature profile for p = 
1.0 x 10 7 g cm -3 , T max = 2.8 x 10 9 K and T = 1.0 x 10 7 K, does not lead to a detonation, 
no matter how large the heated region is. While a runaway occurs and quickly consumes 
all the carbon and oxygen, the re sulting shock wa ve is too weak to form a detonation (see 
also table 1 and the discussion in iKhokhlovi (119911 )). In fact, owing to the lack of gradients 
in induction time, the top hat profile turns into a problem of direct initiation. The heated 
region runs away isochorically, turning the initial conditions into a Sedov problem. The 
strength of the generated shock wave is independent of the size of the heated region, leading 
to the result that for a top-hat temperature profile, a larger heated region (once above a 
minimum critical size) will not make a detonation more probable. In fact, for a while it was 
argued that for this reason the initiation of detonations in degenerate 12 C/ 16 mixtures is 



impossible (e.g. iNomoto et al.lll976l ; iMazurek et al.l 119771 ). This claim, however, was based 



on the assumption of direct initiation. It was shown later that formation of detonations in 
degenerate 12 C/ 16 matter was s till possible via the spontanous mod e of initiation, e.g. from 
a suitable temperature gradient (IBlinnikov fc Khokhlovlll986l . 119871 ). 



In contrast, for the same parameters a linear temperature profile leads to detonation 
at R cr u = 1-2 km (see table [1]). This example explicitly demontrates that there is no 
such thing as a "critical mass" of fuel that, if heated above a certain temperature, would 
necessarily result in a successful detonation. Instead, the details of the temperature profile 
are essenti al for the outcome of t he initiation problem. Therefore, the radii and masses of 
table 1 in lArnett fc Livnei (119941 ) are to be interpreted as lower bounds, constraining the 
conditions for a successful thermonuclear runaway - a necessary, but not sufficient, condition 
for the initiation of a detonation. Indeed, the critical radii for p = 1.0, 3.0 x 10 7 g cm -3 
and T max = 2.0 x 10 9 K obtained here are about an order of magnitude larger. (The value 
To = 2.0 x 10 8 K used in their work is bracketed by two of the choices for the background 
temperature used in this work To = 1.0 X 10 7 K and To = 4.0 x 10 8 K.) 
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3.1. Problem setup and method 



This work uses the flash code (IFryxell et al.l 120001 ) for the numerical experiments. 
flash is an Eulerian code with adaptive mesh refinement (AMR) capabilities that solves the 
compressible reactive Euler equ ations with a directionally s plit implementation of the piece- 
wise parabolic method (PPM) (jColella &: Woodward! Il984l ). The Riemann solver is imple- 
ment ed in a way to handle a general non-polytropic equation of state (EOS) (jColella fc Glaz 
19851 ) . The EOS is appropriate for the compositions, densities, and temperature s encoun- 
tered here (ITimmes fc Arnettl Il999i ; iTimmes fc Swestyl |2000| ; IFryxell et al.l 120001 ) . For the 
nuclear energy release we use an inexpensive 13 species a-chain plus heavy ion nuclear re- 
action ne twork containing 4 He, 12 C, 16 Q, 20 N e, 24 Mg, 28 Si, 32 S 36 Ar, 40 Ca, 4 4 Ti, ^ Cr, 52 Fe, 
and 56 Ni (jTimmes 1999 ; Timmes et al. 2000a ). As suggested in Fryxell et al. ( 1989 ). nuclear 
burning was suppressed in shocks, since the real shock is much thinner in spatial extent com- 
pared to the one spread out over 2-3 zones by PPM. The shock was defined as a region of 
compression (negative velocity divergence) and a significant pressure jump (AP/P > 1/3). 



The time step was chosen to be min(t hydro, 
E int is the specific internal energy [erg g -1 ], e nUl 
rate [erg g _1 s" 1 ] and CFL is the Courant-Friedrichs-Lewy number (jCourant et al. Ill928l ). 



n ) = mm(0.8 CFL, 0.01 E int /e nuc ), where 
is the specific nuc lear energy generatio n 



Limiting the energy added to a computational cell in one time step to 1% of the internal 
energy of that cell is supposed to improve the feedback of the nuclear burning on the hydro- 
dynamics, which are handled in an operator split way. The factor 0.01 is chosen somewhat 
arbitrarily; a smaller factor would be more conservative, but also make the simulations com- 
putationally more expensive. A choice of time step purely based on the CFL number would 
ignore important coupling of the nuclear burning to the hydrodynamics and facilitate the 
spurious initiation of detonations. 

The initial conditions consisted of constant density and composition throughout the 
computational domain. The spatial extent of the computational domain was chosen to 
be 64 x R, which were empirically found to be large enough to cover the locus of oxygen 
ignition for all densities, compositions, and geometries considered. 64 top level blocks with 16 
computational zones each and 7 levels of refinement were consistently used. This particular 
setup means that R was resolved equally in all the runs, whereas the spatial resolution 
(-R/1024) varied in absolute terms from case to case. The simulation time was dependent on 
the particular density, but most simulations were run out to 50,000-100,000 time steps. 

Once a radius R + , for which a detonation ensued, and a radius R~ , for which failure 
was the result, were determined, the critical radius was found by bisection. A trial radius 
Rtry = R+ % R was chosen, and if it led to success (failure) then it replaced the old R + (R~)- 
R + /R~ < 1.1 was the minimal goal, but often the radii were determined to 2 significant 
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figures. This means that each determination of a critical radius involves ~ 10 runs. The 
smallest (largest) R + (R~) is listed as R c ,max (Rcmin) m the tables. 



A code to code comparison with iNiemeyer &: Woosleyl (119971 ). who used the Lagrangian 
hydrodynamics code KEPLER with a 19 nuclide network and comparable zoning, was per- 
formed for their three published low density cases (p < 10 8 g cm -3 ) that have a 50/50 
12 C/ 16 composition. Just as in their work, T max = 3.2 x 10 9 K, T = 4.0 x 10 8 K and a 
linear temperature gradient were selected. In all three cases, the crit ical radii determined 
here agree exquisitely with those found by INiemeyer fc Woosleyl (119971 ) (see table EJ), whose 
Lagrangian calculations were also unresolved. 



3.2. Criteria for success or failure 

The method described above for determining the critical radii relies on a criterion for 
success or failure. The decision was made by inspecting the thermodynamic, hydrodynamic 
and compositional profiles of the initiation simulations at late times. For radii near criticality, 
the evolution of the sub- and supercritical cases initially proceeds in a very similar way. 
The temperature rather uniformly rises in the innermost (originally most pre-heated) region 
where 12 C is depleted first (see figs. |5]and[3]). The rapid almost isochoric burning in the thin 
boundary layer surrounding the peak of the temperature inhomogeneity leads the formation 
of a shock (see figs. H] and [2]), which further steepens via the SWACER mechanism (see 
section [2T2T) . The shock then propagates down the temperature gradient, immediately trailed 
by a narrow region of carbon burning (see fig. [5]), accelerating to detonation speed. 

While the early evolution of the barely sub- and supercritical cases is nearly identical, 
the carbon burning region fails to couple to the leading shock and falls behind (see fig. [5] 
left). Once the advancement of the leading edge of the carbon burning zone had ground to 
a halt the detonation was declared a failure. A success is relatively more difficult to declare, 
as the decoupling of the shock and the carbon burning region may occur well outside the 
originally pre-heated region (see fig. [5] left). No "failures" however were ever observed to 
occur after the successful ignition of oxygen, which due to the longer burning time scale is 
delayed compared to carbon burning (see the space time diagram fig. [6]). Oxygen ignition 
was therefore taken as the criterion for success. At the lowest densities (p = 10 6 g cm -3 ), 
the curious situation arises that the temperature reached after carbon burning falls below 
the oxygen ignition threshold. Detonations propagating steadily for distances exceeding the 
size of the star are observed that fuse the carbon and the carbon burning product neon, but 
fail to fuse oxygen (see fig. |7j). If the carbon burning region remained coupled to the shock 
for for a distance that approached a considerable fraction of the size of the star (~ 2000 km), 
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then the initiation was also considered a "success" . 



3.3. Resolution study 

Ideally, a resolution high enough to guara ntee convergence wo uld have been used. Ex- 



trapolating somewhat freely from the work by iHwang et al.l (120001 ). who perfomed a study 
on numerical resolution of detonation waves using a higher order essentially non-oscillatory 
hydrodynamics scheme, it was expected that convergence may be reached no sooner than at 
a resolution of at least 40 zones per carbon burning length scale. A resolution that high was 
unfortunately prohibitively expensive for the large number of detonation initiation cases in- 
tended for study. The critical sizes presented in this paper were determined with simulations 
that left the region of carbon burning unresolved. To get at least some idea about the impact 
of unresolved scales on the critical sizes, a resolution study was performed. For the resolution 
study somewhat arbitrarily a 50/50 12 C/ 16 O composition with a linear temperature gradient 
at a density of 10 7 g cm -3 , a peak temperature T max = 2.4 x 10 9 K and a background tem- 
perature To = 1.0 x 10 9 K in planar geometry were chosen. For these conditions, the carbon 
burning length scale is a few cm (see the resolved calculation depicted in fig. [8]). The critical 
radii increase with increasing resolution and appear to reach an asymptotic value once the 
carbon burning length scale is resolved (see tabled]). The main thing to take away from this 
study is that the critical radii determined here with unresolved simulations are, while at the 
right order of magnitude, not very accurate. In fact, the critical sizes of a corresponding 
resolved simulation appear to be a factor of two or three higher. 



4. Results 

In this section, the results for the determination of the critical radii (decay constants) 
are presented. To highlight the trends with and dependences on density, temperature, com- 
position, and geometry, they are separated and described in turn. 



4.1. Fuel density 

As expected, the critical sizes for spontaneous initiation of detonations decrease with 
increasing density (see fig. [H]). The dependence on density is quite strong, with a density 
increase of less than an order of magnitude (going from 5.0 x 10 6 g cm~ 3 to 3.0 x 10 7 g cm -3 ), 
leading to a decrease in the critical radii around (and sometimes exceeding) two orders of 
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magnitude (see table HJ). 



Contrary to a statement in iRopke et al.l (l2007bl ). we find that in 12 C/ 16 material 
detonations can form at densities as low as p = 10 6 g cm -3 . However, to get a detonation 
at such low density from gradients smaller than the size of the star, it is necessary to have 
suitably high background temperatures (see tabled]). Unlike at the higher densities, oxygen 
ignition is not required for success. The resulting detonations appear to propagate steadily 
for distances exceeding several thousand km, curiously only consuming the carbon without 
burning the oxygen (see fig. [7]). An initially overdriven detonation that also burns the oxygen 
initially eventually relaxes into a steady propagating detonation that only burns carbon at 
this low of a density. 



4.2. Temperature profile 

4-2.1. Peak and background temperature 

Higher peak temperatures decrease the critical sizes for detonation. However, there is 
a limit to this trend. For a given T and p, there is a T max above which the critical sizes 
appear to asymptote (see tables 0] and [I]). 

Four background temperatures were chosen for the para meter study. They a r e Tp = 
1.0 x 10 7 K (very cold), T = 4.0 x 10 8 K (to compare with iNiemever fc Wooslevl (119971 ). 



cold), T = 1.0 x 10 9 K (hot), and T = 1.5 x 10 9 K (very hot). The effect of the variation 
of the background temperature on the critical sizes is two-fold. 

1. Geometrical. For the linear profile for example, for a given T max the same slope clearly 
corresponds to different radii. 

2. Raised internal energy. A higher background temperature means that the internal 
energy of the fuel is higher, and a weaker shock is needed to ignite it. 

Both effects go the same way; the larger To, the smaller the critical radius. For the lower 
peak temperatures, the dominating effect is geometrical, since the slope of the critical profile 
is independent of To (see table H|). For the higher peak temperatures, the critical slopes 
increase with increasing T . 
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4-2.2. Functional form 

The functional form of the temperature profile can have a very large influence on the 
outcome. The profiles can be separated into two categories: 

1. Sharply peaked. Linear and exponential fall into this category. 

2. Flat topped. Gaussian and higher order Gaussian-like make up this category. 

The critical exponential profile is such that the slope at the origin is very nearly equal (or 
slightly larger) to a critical linear profile for the same T max and T (see figs. [TUl [EH and 
table [5]). This explicitly demonstrates that it is not the amount of mass heated above a 
certain temperature that is the deciding factor but rather the gradient in induction times. 

The critical flat topped profiles behave differently from the critical sharply peaked pro- 
files (see tables |6] and [7j). For small T max , the flat topped profiles more readily lead to 
detonation (see fig. [10]). In fact, for a Gaussian profile, even peak temperatures as low as 
T m ax — 1-6 x 10 9 K (for p = 10 7 g cm -3 ) initiate a detonation at relatively small spatial scales 
of ~ 10 km ( see tab le El), which is far lower than the limit of T max = 1.9 x 10 9 K given by 



Ropke et al.l (l2007bl ). For large T maX} the flat topped profiles require larger spatial scales to 
lead to detonation (see fig. E]). At high temperatures, the flat topped profiles are disfavored 
due to their closer resemblance of a top-hat isochoric blast. The almost instantaneous run- 
away resulting from the high initial temperatures well above the threshold of rapid carbon 
fusion in the flat topped region leads to a synchronized burning of the carbon. The inability 
of such a constant volume explosion in CO matter to directly initiate a detonation is the 
reason why the flat topped profiles need to be "stretched" out until the gradient mechanism 
can operate in the wing of the profile (see fig. [TT]) . The low peak temperature flat topped 
profiles do not share this problem since the initially small temperature gradient in the flat 
topped region is amplified during the longer lasting "smoldering phase" leading up to the 
runaway. 



4.3. Composition 

4.3.1. 12 C/ 12 

The mass fractions of 12 C and 16 O in putative SNe la progenitors is uncertain averaged 
over the star, and allows for even more variation near the surface. To ascertain the depen- 
dence of the detonability on the 12 C to 16 O ratio, we have determined a few critical radii for 
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carbon rich compositions consisting of 70% 12 C, 30% 16 and 60% 12 C, 40% 16 as well as 
oxygen rich compositions consisting of 40% 12 C, 60% 16 and 30% 12 C, 70% 16 in addition 
to the fiducial 50% 12 C and 50% 16 0. As expected, the lower the carbon concentration, the 
larger the critical radii (see fig. [T21 and tables [HI EH HH1 ED- The important and perhaps sur- 
prising result is, that for variations in the carbon fraction within the margin of uncertainty 
as derived from stellar evolution calculations, the critical radii are more strongly dependent 
on the 12 C to 16 ratio than they are on other factors, such as geometry (see section I4~4l) . 



4-3.2. Admixture o/ 4 He 

Critical radii for initiation decrease as the detonability of the fuel is increased. Simi- 
larly, initiation temperatures and densities are much lower if 4 He is added to the fuel (see 
table [12]). The composition of the helium rich fuel considered here consisted of 14% 4 He, 43% 
12 C and 43% 12 by mass, corresponding to a composition of roughly one 4 He nucleus for 
every 12 C nucleus. This represents the most favorable condition for rapid energy release via 
the 12 C(a,7) 16 channel and can be considered somewhat of a "sweet spot" of 4 He admix- 
ture. For this composition, peak temperatures as low as 1.0 x 10 9 K at p = 1 x 10 6 g cm -3 
correspond to critical radii on the order of a few tens of km (see table [T2|). 



4.4. Geometry and curvature 

The differences in critical radii for the spontaneous initiation of unresolved detonations 
from a temperature gradient for spherical and planar geometries is moderate. As expected 
(see section [2T6T) . detonations are more reluctant to emerge in spherical geometry. For the 
cases considered here, the ratio of the critical sizes in spherical to planar geometry is on 
the order of a about two or three (see table [T]). Curvature effects for gradient initiated 
detonations thus have a lesser impact on the critical radii compare d to the direct initi ation 



problem, where the radii are increased by a factor of around 10 (e.g. iHe fc Clavinlll994l ). We 
attribute the reduced importance of curvature in this work to the rather different initiation 
circumstances in the two scenarios. The direct initiation problems starts with a very small 
region of space of extreme over-pressure, a single zone disturbance. Curvature effects at 
early times on the shock of the blast wave are large due to the small radii, i.e. initially 
Rc ^$ Ideti that is, the radius of curvature is smaller than a characteristic length-scale (for 
example the carbon burning length-scale) of a fully developed detonation. In the problem 
of gradient initiated detonations considered for this work, the shock forms further out, at 
distances where R c » ^ et , leading to a weaker influence of geometrical divergence on the 
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critical initiation parameters. 



5. Conclusions 



We have presented ID results of the determination of critical (smallest) spatial scales 
for the gradient initiation of detonations in white dwarf matter. We have quantified and 
tabulated the critical sizes for a range of peak and background temperatures, densities, 
geometry, composition and functional form of the temperature hot-spot. In particular, we 
have shown that the spatial scales required for the initiation of a detonation is crucially tied 
to the functional form of the temperature profile. 

An outcome of the work presented in this paper is tha t the current use of detonation 



conditi o ns in the literature sho uld b e more carefully app lied. iFink et al.l (120071 ). iRopke et al. 



(j2007bl ) , iTownsley et al.l (120071 ) , and IJordan et al.l (120081 ) all rather loosely decide whether or 
not detonation conditions are reached based solely on the peak temperature reached in the 
simulation above a certain density. The sensitivity of the critical radii on the background 
temperature and functional form of the profile shows that caution should be exercised when 
employing such a metric as an authoritative answer to the initiation question. Absolute 
limits derived from linear profiles and cold backgrounds such as the minimum temperature 
or density required for detonation are lowered significantly when taking the possible variation 
of these parameters into account, possibly turning simulations that were declared "failures" 
into "successes" and vice versa. It should be emphasized that choosing any such metric relies 
on the high temperature to occur in unburned fuel, a stipulation which, by the very reactive 
nature of the fuel, should be at least made plausible to exist. 

While this work constitutes an improvement over past efforts, there are several caveats 
to keep in mind when interpreting the results: 



• The simplified 13-species a-chain network is only good within ~ 20% of the energy 
generation rate of a large network. Since the initiation process is highly nonlinear, a 
more realistic network may give different results. 

• The perfectly smooth initial conditions are unrealistic. The effects of noisy initial 
conditions are unknown, but possibly disturb the formation of detonations. 

• The critical radii presented here result from simulations where the carbon burning is 
unresolved. While the resolution study (see section [331) gives some confidence that the 
corresponding resolved critical radii are merely a factor of two or three larger, there is 
no guarantee. 
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• Real detonations are three dimensional. 

In spite of all these uncertainties, the general trends and conclusions, briefly summarized 
below, are likely real and should persist. 

• For the same temperature profile, detonation is more likely at higher density. 

• The dependence on geometry (planar vs spherical) is relatively weak. 

• Composition, especially the presence of 4 He, has a large effect. 

• The background temperature is important (mainly for geometrical reasons) for the 
critical size determination. 

• The functional form of the temperature profile is a vital component of the initiation 

process. 

It should be clear that, instead of relying on tabulations of detonation conditions such 
as the ones listed above, it is clearly preferable to resolve the gradients that initiate the 
detonation in the supernova simulation. There is, however, no hope of resolving the carbon 
burning length scales in a full star simulation in the near future. The dependence of the 
detonation conditions on the functional form of the temperature profile, which is unresolved 
in all extant mult i- dimensional SNe la explosion simulations, warrants that definitive and 
authoritive statements about the failure or success of the initiation of a detonation in such 
simultions should be issued with caution. The hitherto unknown effects of both non-idealised, 
noisy initial conditions (for example non monotic temperature profile or non-zero initial 
velocities) and multi-dimensional instabilities on the critical detonation conditions moreover 
comprise uncertainties that should be elucidated in the future. 
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B523820 to the ASC Alliances Center for Astrophysical Flashes and in part by the National 
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Fig. 1. — Initial temperature profiles considered in this work. 
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Fig. 2. — Shock formation during the early stages of the runaway. The density is 10 7 g 
cm -3 in planar geometry for initial conditions that successfully initiate a detonation. The 
linear temperature disturbance has R = 1.2 km, T max = 2.4 x 10 9 K and T = 1.0 x 10 7 K. 
See also figs. [3j 0] and [5j 
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Fig. 3. — 12 C mass fraction in a space time diagram during the early stages of shock forma- 
tion. Clearly visible is the initial almost synchronous burning of a narrow boundary layer, 
followed by a smeared out reaction front which accelerates to form the detonation. The 
initial conditions are as in fig. [2j 
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Fig. 4. — Pressure during the early stages of shock formation. Visible is the increase in 
propagation velocity of the location of the pressure jump between ~ 60 m and ~ 120 m. 
The initial conditions are as in fig. [2j 
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Fig. 5. — Shown are time series of snapshots of temperature, pressure, mass fraction of 12 C, 
and nuclear energy generation rate at a density of 10 7 g cm -3 in planar geometry for initial 
conditions that fail to lead to a detonation (R = 1.1 km, left) and initial conditions that 
successfully initiate a detonation (R = 1.2 km, right). The linear temperature disturbance 
has T max = 2.4 x 10 9 K and T = 1.0 x 10 7 K. 
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Fig. 6. — Space-time diagram of the mass fraction of 16 O for a succesful detonation in 
spherical geometry for an initial composition of equal parts 12 C and 16 by mass. The 
radius of the initial linear hot-spot is R = 1.6 km and T max = 3.2 x 10 9 K, T = 4.0 x 10 8 K, 
p = 10 7 g cm -3 . The increase in detonation speed is clearly visible as oxygen ignites after 
~ 2 ms. 
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Fig. 7. — Shown are mass fraction of select nuclear species as a function of distance be- 
hind the leading shock. The fuel density is p = 10 6 g cm -3 . The detonation successfully 
propagates consuming the carbon fuel, producing intermediate mass elements (more oxygen, 
magnesium, silicon and sulfur). Their mass fraction appear to reach asymptotic levels far 
behind (~ 100 km) the leadin g ed ge of the d etona tion. Similar results were obtained by 



Imshennik fc Khokhlovl fll984T ) and iKhokhlovl fll989l ). 
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Fig. 8. — Resolved 12 C burning zone trailing the leading shock in a fully developed detonation 
for p = 10 7 g cm -3 and To = 1.0 x 10 9 K. The data are from the run with 15 levels of 
refinement listed in table [3j There is one cross symbol for each computational cell. 
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Fig. 9. — Variation of critical radii with density, all other variables held fixed. Chosen at 
random here T max = 2.0 x 10 9 K, T = 1.0 x 10 9 K and a linear temperature profile in 
spherical geometry. 
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Fig. 10. — Critical temperature profiles for T max = 1.8 x 10 9 K and T = 1.0 x 10 9 K. The 
density is 10 7 g cm -3 . 
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12 



C mass fraction 



Fig. 12. — Variation of critical radii with composition, all other variables held fixed. Chosen 
(at random) here T max = 2.4 x 10 9 K, T = 1.0 x 10 9 K and a linear temperature profile in 
planar geometry. 
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Table 1. Critical Radii for Linear Profile in Planar geometry 



1 \ • 1 

Density 


Tmax 


T 


7") 


7") 


TDSph j TD'plCLTL 


lo s( jSc,l„* ) 


[ 10 6 g cm- 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




[log(K cm^ 1 )] 


1 
1 


Z.O 


1UU.U 


ooU,UUU 


onn nnn 
oyu,uuu 




o.OO 






ir;n a 

10U.U 


OA OOO 


9^ OOO 




A 79 


10 


2.8 


1.0 


1,100 


1,200 


1.9 


6.37 






40.0 


860 


870 


1.8 


6.44 






100.0 


280 


290 


2.0 


6.79 






150.0 


62 


63 


3.0 


7.31 




2.4 


1.0 


1,100 


1,200 


1.9 


6.30 






40.0 


820 


830 


1.8 


6.38 






100.0 


390 


400 


2.5 


6.54 






150.0 


240 


250 


2.4 


6.56 




2.0 


1.0 


6,700 


6,800 


2.6 


5.47 






40.0 


5,400 


5,500 


2.7 


5.46 






100.0 


3,500 


3,600 


2.8 


5.44 






150.0 


1,400 


1,500 


3.1 


5.52 




1.8 


1.0 


32,000 


33,000 




4.73 






40.0 


25,000 


26,000 




4.73 






100.0 


11,000 


12,000 




4.82 






150.0 


3,000 


3,100 




4.99 



Table 2. Code to Code Comparison 



Density 


i?NW 


Rms 


[10 7 g cm" 3 ] 






1 


1.0 - 2.0 km 


1.5 - 1.6 km 


3 


25 - 50 m 


38 - 40 m 


10 


1.0 - 2.0 m 


1.40 - 1.45 m 
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Table 3. Resolution Study 



Density 


T 

' max 


To 




^c,max 


AMR levels 


Resolution 


[10 6 g cm- 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




[cm] 


10 


2.4 


100.0 


390 


400 


7 


39.1 








770 


780 


9 


19.0 








850 


860 


11 


5.2 








1,000 


1,100 


13 


1.7 








1,000 


1,100 


15 


0.4 








1,000 


1,100 


17 


0.1 
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Table 4. Critical Radii for Linear Profile in Spherical Geometry 



Density T max 


To 






! (T -t 0) 

-tt-c max 


[10 6 g cm- 3 ] [10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 


[log(K cm" 1 )] 


5.0 3.2 


1.0 


21,000 


22,000 


3.17 




40.0 


14,000 


15,000 


3.29 




100.0 


4,200 


4,300 


3.71 




150.0 


1,200 


1,300 


4.13 


2.8 


1.0 


20,000 


21,000 


3.13 




40.0 


14,000 


15,000 


3.22 




100.0 


4,100 


4,200 


3.64 




150.0 


1,200 


1,300 


4.02 


2.4 


1.0 


21,000 


22,000 


3.05 




40.0 


13,000 


14,000 


3.17 




100.0 


4,700 


4,800 


3.47 




150.0 


3,100 


3,200 


3.46 


2.0 


1.0 


80,000 


81,000 


2.39 




40.0 


65,000 


66,000 


2.39 




100.0 


40,000 


41,000 


2.39 




150.0 


20,000 


21,000 


2.39 


7.0 3.2 


1.0 


6,900 


7,000 


3.66 




40.0 


4,700 


4,800 


3.77 




100.0 


1,500 


1,650 


4.15 




150.0 


470 


480 


4.55 


2.8 


1.0 


6,800 


6,900 


3.61 




40.0 


4,700 


4,800 


3.70 




100.0 


1,650 


1,700 


4.03 




150.0 


480 


490 


4.43 


2.4 


1.0 


6,800 


6,900 


3.54 




40.0 


4,600 


4,700 


3.63 




100.0 


2,100 


2,200 


3.81 




150.0 


1,300 


1,400 


3.82 


2.0 


1.0 


38,000 


40,000 


2.71 




40.0 


30,000 


32,000 


2.71 




100.0 


19,000 


20,000 


2.71 




150.0 


9,000 


10,000 


2.72 
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Table 4- 


-Oontinupd 






Density T max 






Rc,max 


rp rp 


[10 6 g cm- 3 ] [10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 


[log(K cm" 1 )] 


inn q o 
1U.U 6.1 


l.U 


o onn 
Z,ZUU 


z,zoU 


A 1 £1 

4.10 




,4 a a 
41). U 


1 C AA 

1,5(JU 


l,ol)U 


4. ZD 




1 AA A 

1UU.U 


C O A 


E A A 

59U 


/I CO 

4.5o 




1 C A A 

150.0 


1 DA 

180 


1 AA 

190 


A A/" 1 

4.96 


9 8 

. . . z.o 


l.U 


9 9^0 


9 300 


4 OQ 




40 


1 500 


1 600 


4 1 Q 




1 00 


600 


61 


4.47 




1 (JU.U 


1 80 


1 QO 




2.4 


1.0 


2,200 


2,250 


4.03 




40 


1 400 


1 500 


4.14 




1 DO O 




1 000 


4 1 7 




1 ^0 


^Q0 




4 1 8 


2.0 


1.0 


17,000 


18 000 


3.06 




aci o 


1 4 000 


1 ^ 000 


3 04 




1 00 

1UU.U 


Q 000 


1 000 


3 09 




1 50 


4 500 


4 600 


3 04 


on n o o 

oU.U 6.1 


1 n 
l.U 


A 1 
4 i 


/I Q 


O.OO 




A A A 

40.0 


OA 

39 


/I A 

40 


cr o e 

5.85 




1 AA A 

100.0 


OA 

29 


OA 

30 


5.87 




1 C A A 

150.0 


14 


15 


0.13 


2.8 


1.0 


54 


55 


5.71 




40.0 


46 


47 


5.71 




100.0 


34 


35 


5.72 




150.0 


15 


16 


5.92 


2.4 


1.0 


185 


192 


5.10 




40.0 


150 


160 


5.11 




100.0 


100 


110 


5.12 




150.0 


66 


67 


5.13 


2.0 


1.0 


1,700 


1,800 


4.06 




40.0 


1,400 


1,500 


4.04 




100.0 


900 


1,000 


4.02 




150.0 


480 


490 


4.01 
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Table 5. Critical Decay Constants for Exponential Profile in Planar Geometry 



Density 


-L max 


To 


Rc.min 


Rc,max 


jDexp j rylin 


[10 6 g cm- 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




10 


2.4 


1.0 


890 


900 


0.75 






40.0 


610 


620 


0.75 






100.0 


340 


350 


0.87 






150.0 


200 


210 


0.84 




1.8 


1.0 


28,000 


29,000 


0.88 






40.0 


21,000 


22,000 


0.85 






100.0 


9,900 


10,000 


0.83 






150.0 


2,800 


2,900 


0.93 



Table 6. Critical Decay Constants for Gaussian Profile in Planar Geometry 



Density 


T 

± max 


To 


Rc,min 


Rc,max 


jDgauss 1 jylin 
^c^max! -^c,max 


[10 6 g cm" 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




1 


2.4 


100.0 


140,000 


150,000 




10 


2.4 


1.0 


1,000 


1,100 


0.92 






40.0 


750 


760 


0.91 






100.0 


300 


310 


0.77 






150.0 


110 


120 


0.48 




1.8 


1.0 


9,600 


9,700 


0.29 






40.0 


7,100 


7,200 


0.28 






100.0 


3,200 


3,300 


0.27 






150.0 


1,100 


1,200 


0.39 




1.6 


100.0 


8,800 


8,900 
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Table 7. Critical Decay Constants for glO in Planar Geometry 



Density 




To 


Rc,min 


Rc,max 


r>gW ITtlin 
± *'c 1 max/ ±v c,max 


[10 6 g cm- 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




10 


2.4 


1.0 


1,300 


1,400 


1.17 






40.0 


960 


970 


1.17 






100.0 


250 


260 


0.65 






150.0 


67 


68 


0.27 




1.8 


1.0 


1,500 


1,600 


0.048 






40.0 


1,000 


1,100 


0.042 






100.0 


750 


760 


0.063 






150.0 


570 


580 


0.19 



Table 8. Critical Radii for Linear Profile in Planar Geometry 70% 12 C 30% 16 



Density 
[10 6 g cm" 3 ] 


T 

- 1 max 

[10 9 K] 


To 
[10 7 K] 


Rc.min 

[10 2 cm] 


Rc.max 

[10 2 cm] 


p 70/30 , £,50/50 
-tic.max j -ttc,max 


10 


2.4 


1.0 


240 


250 


0.21 






40.0 


210 


220 


0.26 






100.0 


140 


150 


0.30 






150.0 


94 


95 


0.38 



Table 9. Critical Radii for Linear Profile in Planar Geometry 60% 12 C 40% 16 



Density 




To 


Rc,min 


Rc,max 


R 60/40 , ,-,50/50 
-tic.max j J~tc,max 


[10 6 g cm" 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




10 


2.4 


1.0 


370 


380 


0.32 






40.0 


320 


330 


0.40 






100.0 


220 


230 


0.57 






150.0 


140 


150 


0.60 
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Table 10. Critical Radii for Linear Profile in Planar Geometry 40% 12 C 60% 16 



Density 


T 

max 


To 




Rc,max 


R 40/60 /p 50/50 
ftc,max / *i>c,max 


[10 6 g cm- 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




10 


2.4 


1.0 


7,000 


7,100 


5.92 






40.0 


5,500 


5,600 


6.75 






100.0 


1,800 


1,900 


4.75 






150.0 


520 


530 


2.12 



Table 11. Critical Radii for Linear Profile in Planar Geometry 30% 12 C 70% 16 



Density 


T 

± max 


To 


Rc,min 


Rc,max 


R 30/70 /R 50/50 
rtc^max 1 f*<c,max 


[10 6 g cm" 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 




10 


2.4 


1.0 


50,000 


51,000 


42.50 






40.0 


46,000 


47,000 


56.63 






100.0 


25,000 


26,000 


65.00 






150.0 


2,300 


2,400 


9.60 
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Table 12. Critical Radii for 14% 4 He , 33% 12 C, and 33% 16 



Density 


1 max 


To 




Rc.max 




[10 6 g cm- 3 ] 


[10 9 K] 


[10 7 K] 


[10 2 cm] 


[10 2 cm] 


[log(K cm" 1 )] 


1 


2.0 


1.0 


7,000 


7,500 


5.42 






40.0 


4,500 


5,000 


5.38 




1.6 


1.0 


8,000 


8,500 


5.27 






40.0 


5,000 


5,500 


5.34 




1.2 


1.0 


14,000 


15,000 


4.90 






40.0 


9,500 


10,000 


4.90 




1.0 


1.0 


32,000 


35,000 


4.45 






40.0 


19,000 


20,000 


4.48 


10 


2.0 


1.0 


58 


61 


7.51 






40.0 


45 


50 


7.51 




1.6 


1.0 


100 


110 


7.16 






40.0 


80 


85 


7.15 




1.2 


1.0 


275 


288 


6.62 






40.0 


188 


200 


6.60 




1.0 


1.0 


600 


650 


6.18 






40.0 


350 


380 


6.20 



